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We give a convenient representation for any map that is covariant with respect to an irreducible 
representation of 5(7(2), and use this representation to analyze the evolution of a quantum direc- 
tional reference frame when it is exploited as a resource for performing quantum operations. We 
introduce the moments of a quantum reference frame, which serve as a complete description of 
its properties as a frame, and investigate how many times a quantum directional reference frame 
represented by a spin-j system can be used to perform a certain quantum operation with a given 
probability of success. We provide a considerable generalization of previous results on the degrada- 
tion of a reference frame, from which follows a classification of the dynamics of spin-j system under 
the repeated action of any covariant map with respect to SU{2). 

I. INTRODUCTION 

Quantum operations, such as the application of a unitary rotation on a qubit or a projective measurement in some 
basis, require some form of reference frame. In most descriptions of quantum experiments, this reference frame is 
considered to be classical and therefore can be treated as non-dynamical and used repeatedly without disturbance. 
The situation is quite different for a quantum reference frame - measurements and interactions can cause an unknown 
disturbance on the quantum frame and, consequently, it can reduce the ability of that system to serve as a reference 
frame in subsequent uses [H [H [31 U) |S] . 

Consider, for example, an apparatus implementing a unitary rotation operation on a qubit. For concreteness, we 
consider this qubit to be spin-1/2 system; however, it could equally well be the polarization state of a single photon, a 
two-level atom, etc. A standard apparatus uses a classical directional reference frame to define the axis about which 
the rotation occurs. Now suppose that due to miniaturization, space constraints, or to improve the speed of the 
operation, the system that serves as the reference direction is not accurately described classically, but instead requires 
a quantum treatment using a Hilbert space of finite size. The quantum operation describing the unitary rotation is 
then an operation conditional on the state of the quantum reference frame. Two effects will occur as a result of the 
quantum nature of the reference frame. First, this conditional operation will not be identical to the classical case 
due to the inherent uncertainty in the direction of the quantum reference frame. Second, upon use, the state of the 
quantum reference frame may experience an unknown disturbance, essentially due to entanglement induced between 
the reference frame and system as a result of the interaction. 

It should be noted that the constraints that lead to the requirement that the reference system be treated quantumly 
are ones that are expected to arise in a quantum computing architecture. For example, measuring devices for such 
a system must couple strongly to specific registers (subsystems), and this could imply a very small device. Ideally, a 
measuring apparatus should be well-modeled classically, to reduce noise, however, that may not be compatible with 
the previous requirement. Therefore, it is important to consider how a moderately sized system decoheres when being 
used in this way. This work is directly applicable to the magnetic resonance force microscopy (MRFM) proposal 
of [SI 13 • In that scheme a magnet on the scale of tens of nanometers in length is placed on tip of a nanomechanical 
resonator, which at very low temperatures may allow single-spin measurements in a solid state quantum computer. In 
this regime the assumption of classical behavior of the device leads to a poor approximation. Similar problems would 
occur in situations in which Hamiltonians must be applied to specific subsystems without disturbing neighboring ones. 

We now formally define the scenario that we are considering. Let p*^") be the initial state of the quantum reference 
frame. Let cr**" be the state of a reservoir composed of n ordered subsystems on which we will sequentially perform 
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operations using the same quantum reference frame. Using the first subsystem, we apply a quantum operation x on 
the joint system p^"-* (g) ui. We require that the map x is rotationally invariant, meaning that it is independent of any 
specific direction in space. Formally, a map x is rotationally invar iant if and only if x[R{^){-)R{^)^] = R{n)x{-)R{ny 
for any rotation il e 5*0(3) (or more generally SU{2) if half-integral spins are considered) of the joint system. (Here, 
R is the unitary representation of the rotation group 5*0(3) on the joint system.) Subsequent to the operation, the 
reduced state of quantum reference frame is p^^' = Tr^^ [x(p^'''' "X* Ci)], where Tr^. denotes the partial trace over the 
i*^ subsystem of the reservoir. The first subsystem of the reservoir is then discarded, and the same operation is 
performed on the second subsystem, but using the updated quantum reference frame p^'^\ We repeat these steps 
with the following subsystems of the reservoir, always using the updated quantum reference frame. The state of 
the quantum reference frame after the z*'' step is recursively p^'"*"^-* = Trrjx(p*'*'' ® "'i)]- -^Y discarding the previous 
subsystems of the reservoir at every step, we assume that the reservoir cannot be used to increase our knowledge 
about the quantum reference frame. 

We wish to define a measure of the quality of the quantum reference frame, and to study how it decreases with 
the number of times the reference frame is used. Defining the quality of a quantum reference can be quite arbitrary 
and depends on the interests of the experimentalist and on the task at hand. In this paper, we use the term quality 
function for any function F that is meant to quantify the ability of the reference frame to perform a particular task. 
We will set general conditions on such functions, and analyze how such functions evolve with the number of uses of 
the quantum reference frame. We will also analyze the longevity of the quantum reference frame, that is, the number 
of times it can be used to perform a certain task before its quality falls below a certain threshold. 

As in [3 H] , a quantum directional reference frame could be used to measure whether a series of spin- ^ particles 
are parallel or anti-parallel to some classical arrow. (See Figure [l]) However, there are a myriad of other possible 
cases, including measuring the angular momentum of a series of spin-j' systems with j' > ^, or using the quantum 
reference spin to indicate an axis about which a unitary rotation operation is performed. 

n 
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FIG. 1: The scenario studied by Bartlett et al. The quantum reference spin is used to measure the direction of a series of 
spin- 1 particles in the completely mixed state by means of a projection onto the J — 5 or J 4- | subspaces. 



In this paper, we generalize the results of Refs. [3] H] in three different ways. First, we consider performing 
operations on a reservoir of qudits (quantum systems with d-dimcnsional Hilbert spaces) instead of restricting to the 
case of qubits (c? = 2) as in [5J Second, we consider all possible rotationally-invariant quantum operations on 
the reference frame and the system, while the analyses in [3j |4| were restricted to very particular kind of interaction 
(specifically, interactions describing measurements). Note that this joint quantum operation could be a measurement, 
a conditional unitary, or any other completely-positive map, provided that it is rotationally invariant. Third, we 
generalize the concept of the quality function, and place general conditions on its form and evolution. In [3,4, only 
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a particular kind of quality function was considered - one which was based on the average probability of a correct 
measurement given a known spin direction - while now we are interested in any possible quality function. 

The paper is organized as follows. In Section |ll] we introduce a set of parameters for describing a quantum reference 
frame, which we call moments, and demonstrate that any quality function must depend only on these moments. We 
then present the key result, Theorem[T] which provides a classification of all rotationally-invariant maps. We give some 
recursive equations describing the evolution of the moments (Theorem [2]). In Section III we discuss the longevity of 
a quantum reference frame for a general quality function under the repetitive application of a rotationally-invariant 
map, by analyzing the dependance of the evolution of its moments. Section IV is dedicated to a pair of examples 
intended to illustrate the power of the techniques developed in the earlier sections of the paper. We explore the 
specific cases of measuring spin-1 particles relative to a quantum reference frame, as well as the use of a reference 
frame to perform Pauli operations on qubits. We conclude with some general comments in Section \V\ 



II. MOMENTS AND FIDELITY FUNCTIONS 

A classical reference frame is a local convention for describing the state of a system. A quantum reference frame 
is a quantum system that carries information about a classical reference frame. We define a quantum directional 
reference system as one which is correlated with a classical direction in space represented by the vector n. The 
quantum system that indicates the direction of our reference axis is taken to be a spin-j system, which has dimension 
d — 2j + 1. As in we consider an arbitrary initial state p^^\h) for the spin-j quantum directional reference 
frame. (That is, the quantum reference frame transforms under rotations according to the spin-j representation.) If 
we additionally assume that the initial quantum reference frame depends only on the vector h, we deduce a symmetry 
condition: if Rj{^) is the spin-j unitary representation of the rotation f2 G SU{2) that transforms h to n', then 
Rj{fl)p^j^\h)Rj{fl)^^ = p^^\n'). That is, the set of possible initial states {pf'\n)} is in one-to-one correspondence 
with the set of directions h. 

Because a vector n has an invariance group - the group of rotations about this axis - the symmetry condition 
implies that the state of the quantum directional reference frame is also invariant under rotations about the n-axis. 
This in turn ensures that p^^\h) commutes with the angular momentum operator Jn in the n-direction, and thus 
pj"'' (n) is diagonal in the basis of eigenstates of Jn . 

We now consider how the state of the quantum reference frame is updated through its use in repeated rotationally- 
invariant operations x which act on both the frame and the reservoir. We assume as in [3] that the reservoir initiates 
in a state that is invariant under rotations. (The completely mixed state is one example of such a state, but there exist 
other non-trivial states).^ The assumptions that the initial state of the reservoir and the joint quantum operation 
are spatially invariant imply that the map ^ used to update the state of the quantum directional reference frame, 
p'j'^^^ = iip^^^) = Trr-[x(Pj*^ <^i)]} is also rotationally invariant. We denote the restricted map ^ on the quantum 
directional reference the disturbance map. It describes the back-action on the reference which occurs as a result of 
the interaction. 

Let Rn{9) be the rotation by an angle around the vector h. Suppose that a rotationally-invariant map ^ is applied 
to density matrix pj that is diagonal in a basis given by the eigenvectors of Jn. Because 

Rnm{Pj)Rn\0) = mn{6)p,Rl\e)) = ^(p,), (1) 

then S,{pj) is also diagonal in the basis given by the eigenvectors of Jn. As a result, the evolution of the quantum 
reference frame under the repeated application of the rotationally-invariant map ^ can be described by 2j-f 1 equations. 



Note that, in ^4), they considered the case where the qubits could also be partially or fully polarized, implying that there was some 
initial correlation between the systems to be measured and the quantum reference frame. Generalizing their results to maps x that do 
not preserve angular momentum and to reservoirs that are not polarized in the direction of the quantum reference or only composed of 
spin- 1 particles is a challenging problem. We hope that the framework wc provide here will be useful in extending our results to this 
general situation. 
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one for each of the eigenvalues of pj. A different set of parameters that are equivalent to the eigenvalues is the set of 
moments given by 

{Tr[p,j|] I l^£<2j}. (2) 

Note that only 2j moments are necessary because the sum of the eigenvalues must be one. The use of moments 
instead of eigenvalues will greatly simplify the analysis of the evolution of the quantum reference frame. 

The main motivation for studying the moments of the quantum reference frame is that any quality function F will 
depends only on these moments. Consequently, the behavior of the different moments will determine the behavior 
of the different quality functions. To see why F depends only on the moments given by equation [2] we first remark 
that any reasonable form of F depends only on the state of the quantum reference, but that it must respect the 
relation F{pj) — F{Rj{n)pjRj{n)^^) for all rotations G SU{2) and state pj diagonal in the basis composed of the 
eigenvectors of Jn. In other words, the quality measure should not be biased such that it favors a quantum reference 
frame that is pointed in any particular direction relative to some external frame. All directions must be equally valid. 
Therefore, F does not depend on the direction of n, but only on the eigenvalues or the moments of pj. Note that the 
set of moments with respect to the direction fi can be written as a function of the moments with respect to any other 
direction — this is simply a change of basis. 

Recall that pj'^' as the state of the quantum reference frame after the k*-^ joint operation with a subsystem of the 
reservoir, and p^'^'' as the initial state of the quantum reference frame. By our previous arguments, for a given state 
u^n q£ ^j-^g reservoir and joint quantum operation Xi the moments of pj'^' must be explicit linear functions of the 
moments of . These results lead us to the general recursion relation, 

Tr[pf4]-E^r^[pf-''4], (3) 

where A] are real coefhcients. 

We now present a theorem which provides a classification of the different rotationally-invariant maps in terms of Lie 
algebra generators. The important question of how to obtain noiseless subsystems for quantum channels constructed 
from generators of a Lie algebra has been first studied in [8J, and various proprieties of Lie algebra channels were 
analyzed in [S]. In contrast to those works, we restrict our attention to SU{2). First, some notation. Let Jk for 
k = x^y, z he the angular momentum operators for some arbitrary Cartesian frame. On a spin-j system, define the 
map 

c(P.)^^ E /'^p.^- (4) 

k£{x,y,z} 

for Pj a density matrix of the spin-j system. Let denote the 7?,-fold composition of ^ for n > 0, and define 

Theorem 1 Any map ^ which is invariant with respect to a spin-j irreducible representation of SU{2) has the form 

= (5) 

n=0 

where are real coefficients. 

The proof isprovided in Appendix [Aj 

Theorem [l| allows us an immediate simplification by restricting the number of coefficients A^^"^ of equation (jsj) that 
are required to calculate the evolution of the reference state. The following theorem limits the number of coefficients 
{A[^^) which determine the evolution of the -moments of a spin-j under the action of a general invariant channel, 
and can be very useful to study the behavior of moments with low degree: 
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Theorem 2 If £ is even, then 



i=0 

and if £ is odd, then 



il+l)/2 



Proof. We first define the z-axis so that it is parallel to the vector n which describes our classical directional reference 
frame. Consider any rotationally-invariant map ^. By Theorem [l] we can write £,{pj) — X^fc'Lo 9feC°'^(Pi) where the 
coefficients Qk are real numbers and ^ is given by equation (Q. Therefore, 

2j 

Tr[f(p,) Jf] ^ Tr[p,e( Jf)] = E 9«Tr[p,C°"(4)]- (8) 

ji=0 

To prove our theorem, it is sufficient to show that C°"(>^i) is a polynomial in of degree £, but where all the powers 
have the same parity as £. Define A = j(j + 1) and define the function 

G{£) EE '-{J.jt'jy - Jyjfj.)- (9) 

Using the commutation relations [Ja, Jt] ~ i^a,b,cJc, we can show that 

C{Ji) = aJi-')J. + G{t), (10) 

and 

Gw=G(^-i)j,+c(jr')-^. (11) 

By induction, it is easy to prove using those relations that C{Ji) is a polynomial in Jz of degree £ where all the 
powers have the same parity as £. Using the fact that if (°"{J^) — J2l=o ^sJl for some complex b^, then C°"+^( ) ~ 
C(X]s ^sJl) — ^sC{Jz) and we can prove by induction that C°"(>^f ) is a polynomial in of degree £ where all the 
powers have same the parity as £. This concludes our proof. ■ 



III. LONGEVITY OF A QUANTUM REFERENCE FRAME 

Suppose we have a microscopic device for performing an operation or measurement on a quantum spin using a 
quantum reference frame. Ultimately, after many uses, the quantum reference frame will need to be reinitialized. 
However, we wish to make as many uses of it as possible before performing this reinitialization, without allowing 
the accuracy to fall below some allowed threshold. To define this accuracy we will pick some quality function that 
is suited to our particular purpose. Recall that any quality function F will depend only on the moments of the Ja 
operator. 

As in [3J 13] , we are interested in the scaling, with respect to Hilbert space dimension, of how many times a quantum 
reference frame can be used before the value of its quality function falls below a certain threshold. We refer to 
this property as the longevity of a quantum reference frame. There are three important features of our analysis to 
highlight. First, whereas the work of |31 S| restricts attention to one particular quality function, we will investigate 
the behavior of arbitrary functions satisfying the invariance relation described in Sec. |TT] Because we consider any 
such function, the longevity of the reference frame can be arbitrary. However, we are able to prove some general 
statements about the decay of the moments, and relationships amongst them, and these results can then be used 
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to infer the behavior of any particular quality function. Second, rather than considering only a particular state of 
the reference frame (in [3 S], they considered the state — \j,j)n{j,j\ because it was optimal for the task at 
hand), we consider an arbitrary state. As such, the initial state pj"^ of the quantum reference frame, as specified 
by its moments, can depend on j in a quite arbitrary way. Finally, we consider arbitrary rotationally-invariant joint 
quantum operations x- In [3, 4J, x was chosen to describe a particular measurement that was optimal for some task. 
The resulting disturbance map that they considered is 

where C{Pj) is given by equation Q. This covariant map was also analyzed in [9]. In [3j, it was shown that the 
number of times a reference frame can be used before its quality function (which in their case depends only of the 
first moment) falls below a certain threshold value scales by a factor j^. 

In the following, we prove a general theorem about the longevity of quantum reference frames. We start with 
Theorem[l] which states that any rotationally-invariant disturbance map can be written as ^ = X]n=o 1n(,°^ ■ However, 
as we now show, if the range and values of the coefficients qn can depend on j arbitrarily, it is impossible to make any 
general statements about the longevity. We identify some natural assumptions for the coefficients, which then allow 
us to present a general theorem. 

First, in general, the number of parameters qn describing a map increases with j\ i.e., there are 2j + 1 parameters 
which can be significant to the evolution. As an example, consider the complete depolarization map x{Pj ® — 
2j+i^2j+i ® 5-^1/2- It is straightforward to show that, for this map, all parameters qn for < n < 2j are nonzero. 
However, inspired by an argument in [4 , consider a rotationally-invariant map x{Pj ® f) that conserves angular 
momentum in the fi direction, meaning that x{Jn) = Jh where Jn is the total angular momentum operator in 
the fi direction. The change of angular momentum of the quantum reference frame caused by such a disturbance 
map ^ cannot be higher than the change in angular momentum of the subsystem a, which is in turn bounded by 
the subsystem's dimension d. Given that the map ('(•) = \ 'Ylik=x y z '^k{')Jk cannot lower or increase the angular 
momentum in any direction by more then one unit, the disturbance map's coefficients therefore satisfy = for all 
n ^ d, where d is the dimension of each subsystem a of the reservoir. Thus, if a rotationally-invariant map conserves 
angular momentum, then we can bound the number of non-zero coefficients qn in a way that is independent of j. 

Second, if the parameters qn can depend on j arbitrarily, then we would not be able to conclude anything about 
the longevity of the quantum reference frame. To understand this, consider the example 

j+k 

x{pj(»<T)^ajPj(»(T+{l-aj) ^ Hj+fe(/?j (8)cr)Hj+fc, (13) 

k=b-k\ 

where Hj+k is a projection into the subspace of total angular momentum j + k. This map conserves the total angular 
momentum of the joint system. The dependence of aj on j will determine directly the rate of the decay of the 
moments in function of j. Because the dependence of aj can be arbitrary, then the rate of the decay of the moments 
can also be an arbitrary function of j. We conclude that, in order to make a statement about the longevity of the 
quantum reference frame, we need to assume more than a bound on the number of the (/„ parameters or that the 
rotationally-invariant map x conserves angular momentum in the direction of h. For the following theorem, we assume 
that each coefficient qn converges to a constant when j — > cxd. We also need to make an assumption about the rate of 
convergence: suppose that each g„ can be written as a quotient of two polynomial in j, such that the degree of the 
denominator is at least the degree of the numerator (i.e., qn ^ 0(1)). Finally, we assume that the state p^^^ has the 
propriety that Tr[pj°-' j|] = 0{j^). A sufficient condition for this to be true is that there exists a /3 > independent 
of j such that j)fi > P- (We note that these assumptions are motivated by examples, which we explore in 

the following section.) 

Theorem 3 (Longevity) Consider a quantum reference frame, realised as a spin-j system with initial state pf'\ 
which is used for performing a rotationally-invariant joint operation x- If this operation induces a disturbance map 
^ = X^nLo ^^'^^ satisfies the following assumptions: 
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1. there exists some Umax such that = for all n ^ rimax , 

2. qn ^ 0(1) , 

3. Tr[pfhi]^0{f), 

then the number of times that such a quantum reference frame can be used before its moment falls below a certain 
threshold value scales as p. 

Proof. Noting that the Jk operators are self-adjoint and using the cychc propriety of the trace, the map C, has the 
property that 

Tr[C(p)4]=Tr[pC(4)]. (14) 

Therefore the moments of angular momentum in the h direction after the map has been applied to the reference 
state can be expressed as 

Tr[pf+^)j|] = X:9nTr[pf C°"(4)] (15) 

Using the commutation relations, the factor C°^{JfJ can be expanded as a polynomial in Jn of degree £. Define 
^ — iU + 1)- Using a proof by induction on equation (10 1 and observing that J2ke{x y z} JkJnJk — ~ ^)Jn, it is 



easy to show that the coefficient of the leading term will be A^p — X]n=o°'[9n(l ~ '^(x))] = 1 ^ ^{\) where we used 
the normalization condition 9" = Ij assumption that each g„ is 0(1) and the fact that q^ = for n > 2j. The 
coefficients of the non- vanishing lower terms will be = 0(1) for i < i. 

This reasoning about the constants A^,*^ is sufficient to characterize the rate of change of the moments with repeat 
application of the map. To demonstrate this, let ^ = 1. We want to find the minimum value of t such that 
Tt[pP Jn] < c, for some constant c. Using equation tl\, we need to solve 



c = Tr[p^."^J^](l-0(-) ) . (16) 



Therefore, 

tc^O {f) . (17) 

To generalize to higher moment, observe that Tr[p'°'' J^] ^ 0(j*) for i < £. Recall that we assumed that Tr[p^°'' j|] = 
0{f). Using equation the facts that A'-p = 1 - 0(i) and A^ = 0(1) for i < £, we can extend our result to 
higher odd moments by using strong induction. In other word, we can show that the minimum value of t such that 
Tr[pf 4] < c for proper chosen initial state Pj'^' is O(j^). For even moments, a similar approach using equation (joj) 
can be used to obtain an identical conclusion. ■ 

Note that in the case were the assumptions of Theorem [3] fail, but there is a specific dependance on j of the map x 
and the state p^^^ , equations ^ and (j? I may still be useful tools to study the longevity on quantum reference frame. 



IV. EXAMPLES 

A. Measuring Spin-1 Systems 

As an example, consider the measurement of spin systems relative to a directional reference frame. This problem was 
investigated in [3J[3] for the case of spin-1/2 systems. Specifically, suppose that the reservoir consists of spin-s systems 
which are initially in the completely mixed state. The quantum reference direction is a spin-j system with j > s, 
aligned with some classical direction n. The goal is to measure each spin to determine its component /i along the vector 
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n by performing a joint measurement on both the system from the reservoir and the reference system. The optimal 
rotationally-invariant joint operation for this task [TU] is a POVM given by the projectors {Ilj+^l/j, S {— s, . . . , s}} 
where Hj+p corresponds to a projection onto the subspace where the total angular momentum of the reference spin 
with the measured spin is j + fi. Define the fidelity Fs as the probability that, when the reference is used to measure 
a system in a known state |s, ii)n by means of the above POVM on the joint system the result it returns is correct, 
assuming that all values of fi are equiprobable. 

Consider briefly the case of s = 1/2, as in [3, 4 . The fidelity function can be expressed in terms of the first moment 
of the reference frame after being used k times, as 



5 + 5-nTrb<«Al 



(18) 



Using Theorem [2] and a simple calculation to evaluate A'"^\ as defined in equation we find the fidelity in terms of 
the original value of the first moment to be 



i 2^ 2j + l V (2j + l)2 



(19) 



Note that the fidelity function in this case depends only on the first moment. 

We now generalize to the case where we wish to measure the angular momentum of a spin-1 particle along some 
direction using a quantum direction reference frame. (Larger spins can be handled by similar means.) First, we 
determine an expression for the fidelity in terms of the moments: 



J^n,+^(pf ®|A*)(/i|) 



where fj, G {—1, 0, 1}. Consider the disturbance map f as described in Section [ll} 



(20) 



(21) 



Explicit expressions for the reduced operators on the reference system (/i"|nj+^/ can be found from the Clebsch- 
Gordon coefficients. The expression for the fidelity in terms of j and the moments is 



_ 1 , [(2j + 1)^-2] (fc) 1 (fc) 2, 

~6 + 6j(j+i)(2,+i)^'^[^^- ''^^m^)^'^'^ 



(22) 



Note that the fidelity in this case depends on the second moment as well as the first. Using Theorem |2] we have: 



TrK(pfVf] = A^> + A^^>TrW'Jt] 



(2) 



,(2)^^r (fe) .21 



(23) 
(24) 



Substituting some particular values of pj (i.e., pj — \ j, m)n{j, m\ for m — j, j — 1) and solving the system of equations, 
we obtain: 



(1) 



3j2(j + l)2 - 3f 

(2) _ 2 (8/ + 16 f - 8j - 3) ^ 4 5 
3j(j + l)(2j + l)2 3 3f 





(2) 



0|i 



4f + 12 f ^ 3j* - 26 f + f + 16j + 6 



(25) 
(26) 
(27) 
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(The series apply as j — > cx).) In terms of these constants, the fidehty evolves with k, the number of applications of 
the map corresponding to the measurement of the z projection of the spin-1 system, as: 



ik)_l [(2j + 1)^ - 2] 
'6 6j(j + l)(2J+l) 



1 



2j(j + 1) 



A, 



(2) 



1 - 



l-A 



(2) 



(28) 




FIG. 2: A plot of the longevity, as defined in Section [III[ against the reference system size j for the scenario described in 
Example A (points converging to the solid line of gradient 2), and Example B, Method 1 (points converging to the solid line 
of gradient 1), Method 2 (circles converging to the dashed line of gradient 1), and Method 3 (grey circles converging to grey 
dot-dashed line of gradient 2). The points are found numerically. The gradients of the lines in A and B(3) are 2, indicating 
longevity scaling as 0{j^), and the gradients of the lines in B(l) and B(2) are 1, indicating a scaling of 0{j). 



Because the assumptions given in Section [III| are satisfied, the longevity of the quantum reference frame must scale 
as 0{p). This result is easily verified numerically. See Fig. [2] 



B. Implementing a Pauli Operator on a Qubit 




The quantum reference frame can also be used to define an axis for the purpose of implementing some desired 
unitary gate. Having only one axis, we are restricted to the set of gates corresponding to rotations about this axis. 
Suppose, for example, we want to implement a Pauli Z operation on a qubit. 



(29) 



using a quantum reference direction pj in the form of a spin-j system to define the z-axis. To do this, we must imple- 
ment a rotationally-invariant operation on the combined reference and system. Unlike the measurement considered 
in previous examples, in this case we wish to implement a unitary operation on the system that is conditional on 
the state of the quantum reference direction. In the language of quantum information, we want to perform a unitary 
rotation about a direction programmed by the state of the reference frame. However, that because the size of the 
reference frame is bounded and the number of programmed operations (possible directions) is continuous, then the 
impossibility of a "programmable quantum gate" [11_ ensures that an ideal such conditional unitary cannot exist. It 
is possible to approximate it, as we now show. 
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1. Method 1 



Consider the following operation, to be performed on the joint state of a quantum reference frame pj and a single 
qubit a. First, define an covariant measurement on the spin-j reference frame described by the POVM 

{A(J7) = (2j + l)i^(^])|e>(e|i?(^7)^ 17gS'C/(2)}, (30) 

where |e) is a normalized state of a spin-j system. The measurement effects of this POVM satisfy the normalization 
condition 

' d^lnm) = Ij , (31) 

where Ij is the identity on the spin-j Hilbert space and dfiQ represents the Haar measure over SU (2) . This measure- 
ment is performed on the quantum reference frame state pj, and then conditional upon obtaining the outcome O, 
the operation Z{il) = R{il)ZR~^{U) is applied to the system. If we subsequently discard the information about the 
measurement result A, then the effective joint map 

x{pj(»a)= [ dfinVH^) ® z{n){pj (g) (T)^/A{n) ® z{n)^ , (32) 

is clearly invariant. (Although the covariant measurement makes use of an external reference frame, the resulting 
map X is independent of this choice of frame.) The net operation on the system a is given by the map 

T{a)= I dnnT^rj[A{n)pj]Z{n)aZ-\n) . (33) 
Jn 

Note that this operation r is also rotationally invariant; although the measurement and subsequent unitary appear 
to require an external spatial reference frame, the net operation is invariant under changes of this frame. Also, 
because the operation is constructed explicitly from a POVM measurement and a unitary operation conditional on 
this classical result, followed by tracing out the state of the reference frame, it is necessarily completely positive and 
trace preserving (CPTP). (This fact is also clear by observation: the term Tr^ [A(ri)pj] is a normalized probability 
distribution weighting possible unitaries Z{fl), and thus the map r is a valid unital CPTP map.) 

This expression for t explicitly gives a Kraus decomposition of this map (albeit with a continuous number of Kraus 
operators), T{a) — J^dpiiE{n)aE{ny , where 



E{n) = ^Trj[A{n)p,]Z{n), (34) 

are Kraus operators satisfying diiQE{^)^ E{iV) = I. The ability of this operation to approximate the Z operation 
on the system is defined using the gate fidelity [121 ESI Ell El] given by 



F.„.(^..).i*^»^ (35, 

where d (the dimension of the system on which the gate is applied) is 2 in our case. 

Suppose that quantum directional reference frame pj is aligned with the z-axis. Then pj is a mixture of states 
\j,m)z for —j ^ TO ^ J. To simplify our calculation, assume that pj = [jjUi) z{j^m\ (i.e., we consider only a pure 
state, but at the end, we can generalize for any mixed state). The state |e) that defines the POVM could be any state, 
but for simplicity we consider only the case where |e) — \j,j)z- In this case, the rotationally-invariant operation on 
the combined reference and system is 

r(a,p,) = (2j + 1) / d^ln\z{j,m\R{n)\jJ),\''z{n)aZ-\n). (36) 
Jn 

The Kraus operators are then given by 

Eifl) = ^{2j + l)Uj,m\R{n)\j,j)JZ{n) . (37) 
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We can parameterize 

\ — e sm e ^ cos / 
where < a,ip,(l) < 2tt and < 6* < ^. The rotation _R(ri) can be rewritten using Euler angles: 

Rin) = e'"i?.(-0 - ^)Ry{-2e)R,i^ - 0), (39) 
where Rk{P) is a clockwise rotation of angle /3 around the /c-axis. It follows that 

Tr[R{n)ZR\n)Z] = Tr[Ry{-2e)ZRy{29)Z] = 2 cos 26*. (40) 

From [11], we have that 



,{j,m\R{n)\j,j), = e~'"^"''^'>'+'l'^-'^^'^-'l''>J(^^'^^{cosey+^ (41) 



Therefore, 



\Ti-[E{n)^ Z]\^ = {2j + l)( .'^^ Vcos0)2J+2"(sin 0)2^-2™ |Tr[i?j,(-20)Zi?j,(20)Z] I' 

= 4(2j + l)f Vcos6l)2J+2™(sin6')2J-2m(^Qg26»)2, (42) 

where, in the first step, we used the cyclic propriety of the trace. The Haar measure in these coordinates is dil.{U{2)) = 
(2tt)-^2 sine cos ededad(l>d^.'^ It follows that 

p (7 ^ 1 , 2(j + l + 2m^) 
^--(^'") = 3 + 3(, + l)(2, + 3) 

which is a function of the second moment of the z projection of the reference frame. The result easily generalizes to 
the case where pj is a mixture of pure states of the form \j, ni) z'. 

tgate{^,T)- ^+ 3(j + l)(2j + 3) • ^^^> 

This example is a case where the fidelity evolves with the expectation value of the second moment, despite the fact 
that the reservoir is composed of spin-i particles. Also note that in this particular case the fidelity does not depend 
on any odd moments because the operation depends only on the axis defined by n, but not on the direction along 
this axis. 

We now consider how the quantum reference frame degrades with repeated use in performing this operation. 
We assume that the qubit systems (the reservoir) on which we apply the approximate phase gate are all in the 
completely mixed state With each application, the reference frame evolves according to the invariant map 

C(Pj) = Trsx(pj (g) ), where Trc, is the partial trace of the qubit system on which we apply the approximate phase 
gate. To calculate how the second moment evolves as a function of the number of times the quantum reference frame 
has been used to perform the approximate phase gate, we can use Theorem [2j 

TrK(p,-) Jf] = + 4''Tr[p, J|]. (45) 



^ We take the integral over the Haar measure for f/(2) instead of 5(7(2) because this simplifies the notation, however, it should be noted 
that in this case the integral over a (the global phase) will not alter the result, so we are free to do this. 
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(2) (2) 

To find the values of the coefficients AJ, and A2 , we consider two possible initial states of the quantum reference 
frame. First, suppose that pj ~ 2jTT^-" ^^i^^ yields (,{pj) = 2j+i -^j- This evolution gives 

ZO+l) ^42)^42)^0+1)^ (46) 
Second, using pj — \j, j) z{j, j\, we obtain a second equation: 



j(l+j(3 + J+2j2)) _ (2) , .(2) ,2 



(l+.7)(3 + 2j) 



= A}:'+A',^>j\ (47) 



Solving those equations, we obtain 



and 

^ (2j- + 3)(j + l) J ^ ^ 

From the above equations, we can deduce that the longevity of the quantum reference frame is 0{j) by noting that 
the fidelity after k repetitions of the gate is 

k 



3 3(j- + l)(2j- + 3) ■ ■ " i_42) 



("2) ("2) 

From the equations (48) and (49) we have that Aq ' goes as 0(j) and A\ ' goes as 0(1). Looking at equation (50) we 
see that fidelity must go as 0(j) and therefore longevity goes as 0{j). This is also seen in the numerical work shown 
in Figure[2j This result appears to be in contradiction with Theorem|3] However, this contradiction is resolved by the 
fact that one of the assumptions of the theorem is not fulfilled. Indeed, even if the map x is rotationally invariant, it 
does not conserve the total angular momentum. In particular, note that (j, TO|z^j(|_7, j|)|j, m) > for all m in 
the range — j, . . . , j. Because C°^"'(|ji i)2(j, j|) is the only map of the form {\j ^ j) z{j ^ j\) for < n < 2 j + 1 that has 
support in the state |j, ^j)z{j, then q2j > 0. Therefore, there is no bound rimax independent of j such qn = for 

all n > TlyyiaX- 



2. Method 2 



We now investigate an alternate method for approximating a Pauli Z operation using a quantum reference frame. 
Consider the filtering operation [T7] 



r = (2j + i) / d^inRJm®Rl/2m \j,j)z{j,j\^z R,{n)-^^R,/2{n)-\ 



Observing that in spin notation 



^ ~ I2' 2)2(2' 2I I2' 2)2(2' 2I 



(51) 



(52) 



we can replace the operator Z in equation (51) with this expression and using the Clebsch-Gordon coefficients we 



can express the part of equation (51) that is in the square brackets as a joint system with total angular momentum 
having support on the j — and j + | subspaces. F is spatially-covariant and we can use Schur's lemma to note that 
it will be block-diagonal in the irreps j + 5 and J — |. Thus, we can rewrite it as 

F = (2j + 1) / d^^nRJ+l/2in) [|j + i, j + + 1, j + i| - + i, j - + l,j - 1 |j i?^-+i/2(r!)-i 

Jn J 

-{2j + l) I d^i^R,^,/2m\^^\J - ^ \)z{j ^ \,j ~ \\\R,-,/2{nr^ (53) 
Jn ^ J 

= (2J + 1) (2^2(1 - 2^1) W2 - 2121^1^.-1/2)' (54) 
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which simplifies to 



r-(2^W2-Vi/2)- (55) 



We can define 

Xj{pj®cj) = T{pj®(j)T\ (56) 

where Ilj+fe is a projection into the subspace of total angular momentum j + k. However, this map is not trace 
preserving, nor does it conserve angular momentum. As in method 1, let pj — |j, m)^ (j, m|. In that case, we have 
Kraus operators of the form 

Era' = z{j,m'\T\j,m),_ (57) 
and, using the appropriate Clebsch Gordon coefficients, we arrive at an expression for fidelity of the form 

1 2 / 1 ^ ^ 



^---3 + iljTTj 

which goes to 1 in the limit j ^ oo, for the case m — j. For a general state pj (which must be diagonal in the basis 
"i)^ I -j^mi^ j}), we obtain 



^^ffate = 3 + 3 ) Tr[p,J,^], (59) 



By Theorem [2] even though the map is not trace preserving, we have: 



and, by the same method as in the previous subsection, we have for this map 

.(2) _ j ^1 1 I 1 

U + 1) j P ' ^ \ f 



TrK(p,) J|] - 4'^ + 4"Tr[p, J,2], (60) 



< = TTTn =1-7 + ^ + 0^)' (61) 



(j + 1) V j(j + i)/ J f w ' 



where the series expansions are appropriate in the limit j — > oo. 

Again we find the scaling of longevity with j to be 0{j). (See Figure [2j) However, the longevity is consistently 
higher than by using Method 1. 



3. Method 3 

We might ask about performing an operation similar to the previous case, but in which trace is preserved on the 
map. Such a map could be realized through a coupling between the spins of the form: 

H^w (ji^) jO) + J(^) + Ji') Ji^)) = f - I + 1) + ^) 1 (63) 

(-) 

where w is some constant, the operator Jx = ^ is the operator on a spin-half system and the operator J^^ — 
(j + i) + I) n^+i + {j — I) (i + ^) is the total angular momentum operator on the joint system of the 

reference and the spin-half system. Therefore, the unitary associated with this Hamiltonian is J7 = e*^^'^ 
Choosing t — ^(^2j+i) gi^^s the appropriate phase shift of the spin-half system relative to the spin-j reference system. 
It could also be phrased in the language of the previous subsection if — up to a global phase — we set 



r = u = (/j+i/2 - ^-1/2)- 



(64) 
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The quality function for this map is 



F, 



gate 



1 2 

3 ^ 3 V2j + 1 



Tr[p,j2], 



(65) 



We have for this map 



which, when Tr[pj J^] = j^, also tends to 1 in the limit j 



A. 



:2) 



1 



^('^ - 1 - 

' (2j + l)^ 



(2j + 1)2 
12 



-> oo. 

= 1 - 
= 1 - 



1 
3 



o 



J 



1 



(66) 
(67) 



where, again, the series expansions are appropriate in the limit j oo. Observe that in this case these parameters 
lack any terms in i . 

Now we find the scaling of longevity with j to be 0{j^), as depicted in Figure [2] Figure [3] compares the fidelities 
for the three methods for constant j. These last three examples demonstrate the importance of choosing the gate 
carefully to match the objectives of the given task. 




FIG. 3: A plot of the fidelity with number of repetitions, n, for j = 8 for the three methods, B(l) (dot-dashed line), B(2) 
(dashed line), and B(3) (solid line). This behavior of this value of j is representative. 



CONCLUSION 



This study has explored the physically relevant question of how a quantum directional reference frame state evolves 
under the action of maps invariant with respect to rotations in space. We generalize the concept of quality of a 
quantum directional reference frame introduced by [3] . We argue the quality of a reference frame must be represented 
by a function that depends only on the eigenvalues of the quantum reference frame or an equivalent set of parameters 
called the moments. We give recursive equations (Theorem [2]) for how the moments evolve with the number of uses 
of the quantum reference frame. We derive sufficient some conditions (Theorem |3]) for the longevity of a quantum 
reference frame to scale by a factor proportional to square the dimension of the quantum reference frame. Finally, we 
applied our results to different examples such as the use of a quantum directional reference frame to measure a spin-1 
particle or to implement an Pauli operator on a qubit. The tools that we developed can be use to compare different 
methods to perform some operation using a quantum reference frame as we showed in our last example. 
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To prove Theorem [2] we use Theorem [T] which impHes that maps invariant with respect to SU{2) can be written 
in a polynomial form as a function of the Lie algebra generators. It would be interesting to investigate if a similar 
theorem could be applied to other Lie groups. Also, we assume in this analysis that the state of the reservoir is 
invariant under rotations. However, in the most general physical situation, the reservoir can be polarized. Building 
on the work done in it would be interesting to investigate how our theorem can be generalized to the case where 
the above symmetry is broken. 
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APPENDIX A: PROOF OF THEOREM 1 



Here, we provide a proof of Theorem [T] inspired by ideas from [HI [19]. First, because the group 5*0(3) of space 
rotations and the group SU(2) of unitary transformation are locally isomorphic up to a phase, then the spatial 
covariance condition can be replaced by covariance with respect to SU{2). 

Let pj be a density operator on the Hilbert space of a spin-j' system. Consider the map 



C(Pi) = jiJxPJx + JyPJy + JzPJz), (Al) 



where A — j{j + 1), and J^, Jy and are the angular momentum operators in the x, y and z directions for some 
arbitrary Cartesian frame. First, we notice that 

ai,)^\jzJl'h^\j'^i,, (A2) 

i—l 

where J is the total angular momentum operator. 

Lemma 1 The map C, is covariant with respect to the spin-j irreducible representation Rj of SU(2). 

Proof. If Rj is the d = 2j + 1 dimensional irreducible representation (irrep) of the group SU(2) generated by 
{Jx, Jy, Jz}, then it can be expressed in the form 

Rj{n) = e-'^^J'^e-^'^^-^ye-'^'-^^ (A3) 

for some real parameters 9i, 02 and ^3 (which are strongly related to the Eulers angles). By the symmetry of C with 
respect to the y and z-axes, if we prove that 

RM'\{R,{e)pjRM~^)Rz{e) = c(p,), V e e [o,27r), (A4) 



for any rotations around the z-axis only, then (A4) with any rotations around the y-axis follows immediately. From 



(A3), we deduce that (A4) is satisfied for any rotations. Such z rotations will map 

Jx cos 6Jx + sin 9Jy (A5) 

Jy —f — sinOJx + cosOJy, (A6) 

and it can be readily computed that ^ is invariant with respect to z rotations. Therefore, 

Rj{n)~'\{Rj{n)pjRj{n)-'^)Rj{n) ^ C{pj) yn£su(2). (A7) 
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The map C is therefore invariant. ■ 

Because a composition of invariant maps is also invariant, then any map of the form 

2i 

^(p) = '?0P + E9fcC°'(p), (A8) 

where the values qi are real numbers and C"*^ = ^ ° C ° • • • ° C; is also invariant with respect to SU{2). It remains to 
show that any invariant map with respect to SU{2) can be written on the form (A8). 

To prove that every invariant map can be written as above, we use the Liouville representation of a superoperator 
[20] , Upon representing a, d x d density matrix p into a d"^ long column vector \p)) by stacking the columns of the 
density matrix, the action of any superoperator ^ can be represented as a x d^ matrix /C(^), such that 

m))^m\p))- (A9) 

This representation is necessarily basis dependent. If a given process has Kraus operators {E^}, the Liouville repre- 
sentation takes the form 



k 

where * represents the complex conjugate with respect to the chosen basis. 

Lemma 2 The Liouville representation of any map that is invariant with respect to SU (2) has the form 

m) - (e""-"" ® 1) ^ Ckllkie^-'y ® 1), (All) 

fe=0 

where G C and Hk is the projector into the 2fc + 1 dimensional subspace of total angular momentum k. 
Proof. The condition on a map to be invariant can thus be expressed as 

o Rj{n))jc{0 = JC{0{R*{n) (g) Rj{n)), \/ne su{2), (A12) 

so the condition on the operators in this new representation is that /C(^) must commute with R* ® Rj. If the group 
SU{2) is represented in the {\j,m)z} basis (i.e. the eigenstate of J^), we have the relation R*{^) — e^"'^''i?(ri)e"'^'' , 
which implies 

{Rj{n) (g) i?j(f]))(e"^'' (g) I)IC{0{e-'''-^y (g) J) = (e"''" (g) /)/C(e)(e-"^'' g) I){Rj{n) g) Rj{n)) Vne SU{2). (A13) 

We note that {Rj{il) (g) Rj{il)) is the collective representation of SU{2) on two spin-j systems. A Clebsch-Gordon 
decomposition gives us the irrep of the group of all collective rotations G {Rj g) Rj ) which take the form 

2i 

g{Rjg)Rj){n)~^Rk{n), ynesu{2), (ai4) 

fc=0 

where "~" denotes "unitarily equivalent" and Rk is the spin-fc irrep of SU{2) which has multiplicity one. 

Because we require /C(^) to commute with Rj{^) <g Rj(S^) for any e SU{2), it must then commute with all irreps 
of G{Rj g) Rj). By Schur's lemma, we have that 

2j 2] 

{e'^'y g> l)/C(^)(e-"^« (g 1) ~ 0Cfe/fc = Jl^^nfe, (A15) 

k=0 k=0 

where Ck € C, Ik is the 2fc + 1 dimensional identity operator, and Hk is the projector into the 2k + I dimensional 
subspace of total angular momentum k. m 

There are 2j + 1 independent projectors forming /C(^). To characterize every possible invariant mapping, we will 
thus also require 2j + 1 independent operators. 
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Lemma 3 The matrices J* ® Ji)^ forO are linearly independent and form a complete basis to represent 

any invariant map of the form (All). 



Proof. If we rewrite equation (A8l in the Liouville representation, the map becomes 

fc 



k=0 



\ie{x,y,z} 



(A16) 



The hermitian matrices {J^i Jt ® Ji)^ are diagonal in the same basis, and the eigenvalues of J* ® Ji)^ are vf , 
where vi is an eigenvalue of ^ ■ J* ® Ji . 

To find all of the eigenvalues {vi}, expand the total angular momentum J'^ — 2- {Ji® I + 1 ® Ji)^ to get a new 
expression for Ji ^ Jf- 

^J,®J, = '^(^{J,® I + I ® J,f - ^Jf ® I - I ®^jf) 



I^J"^) 



-(^2_2j(j + l)/) 



(A17) 



With the relation J't ®Ji^ -e""-^' 
have the same eigenvalues up to a negative sign 
From equation (A17l, we find that 



I /, we thus have that J* ® Ji and Jt ® Ji will 



vi = --{l{l + l)-2]{j + l)), 



(A18) 



where l{l + 1) is the eigenvalue resulting from the vector addition of the two spin-J systems, which implies that I 
can range from to 2j and has multiplicity 21 + \ . There is 2j + 1 different values of so there is 2j + 1 different 
eigenvalues. By the fundamental theorem of algebra, this implies that there exist no polynomial of degree 2] that has 
the eigenvalues {vi} as roots. This impHes the matrices J* (8) Ji)'^ for < fc ^ 2j are linearly independent. By 
counting the number free parameters, we proved that the operators J* ®) Ji)'' form a complete basis to represent 
any map represented by equation (All I. ■ 



Finally, to show that the equation (A8l is a representation of all possible invariant maps on a spin-j system, we 



need to show that the qi's are real. First, note that 



(A19) 



where J± — Jx ± iJy Also, from equation (Al9l, the only contribution to n{j, j)nn{j, j\)\j, —j)h is from 

q2jC^'{\j,j)nn{3,j\)- Actually, 



n{j,-j\^{\j,j)nn{j,j\)\j,-j)n = q2j h{j , - j\C'^' i\j J) fin{j , j\)\j , -j) h ■ 



(A20) 



Because ^{pj) must be positive, then q2j is a non-negative real number. The contributions to n{j,—j + 
m{\j,j)nn{j,j\)\j:-j + l)n are from q2jC'^H\j, j)nn{jj\) and q2j-iC^^~'^{\j, j)nn{j, Since ^{pj) is positive and 
q2j is real, then 92^-1 must also be real. Note that q2j-i could be negative. By induction, it is easy to show that all 
the coefficients qi must be real. 
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